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Abstract 

We study the mean field dilute model of a ferromagnet. We find 
and prove an expression for the free energy density at high tempera- 
ture, and at temperature zero. We find the critical line of the model, 
separating the phase with zero magnetization from the phase with 
symmetry breaking. We also compute exactly the entropy at temper- 
ature zero, which is strictly positive. The physical behavior at tem- 
perature zero is very interesting and related to infinite dimensional 
percolation, and suggests possible behaviors at generic low tempera- 
tures. Lastly, we provide a complete solution for the annealed model. 
Our results hold both for the Poisson and the Bernoulli versions of 
the model. 
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1 Introduction 



The study of the Curie- Weiss mean field ferromagnetic model can be per- 
formed by very simple convexity methods [21 HD] • The main ingredient of 
such approaches is the quadratic dependence of the Hamiltonian on the 
magnetization, which makes the free energy convex in such a quantity. A 
very strong analogy allows to reproduce the same physical approach when 
the ferromagnet interaction is replaced by one with Gaussian distribution 
[21 HO] • In this case the convexity arguments apply through the quadratic 
dependence of the covariance of the Hamiltonian on the main physical 
quantity for these models: the overlap between two configurations. De- 
spite the initial belief that the general approach reviewed in [51 [TU] was 
possible only because of the special nature of Gaussian interactions, it 
turned out that the mentioned analogy extends to dilute mean field spin 
glasses [3 [21 13] , and independetly of the (symmetric) distribution of the in- 
teractions. Surprisingly, the case of dilute mean field ferromagnets has not 
been framed so far in the context of these methods, reviewed in [^ITU1[5]. 
This is not the only reason why the mean field dilute ferromagnet is an 
interesting model, as we will see. It is also noticeable that the mean field 
dilute ferromagnet has not been studied so much. Most of the studies in 
the physical literature are not fully rigorous and only graze the dilute fer- 
romagnet within works about different models, and often within a general 
study of models on various kinds of networks, which gained recently a large 
attention (see for instance [5l[14l[TT] and references therein). An exception 
is [7] , which is a rigorous quite general study about reconstruction for mod- 
els on random graphs, but also analyzes the case of an Ising ferromagnet. 
Our methods and purposes are different from those of [7], where results 
of quite general nature are present about the existence of limits. A clear 
explanation of the physical behavior of the dilute mean field ferromagnet 
is still missing. 

We started our work studying the annealed version of the mean field 
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dilute ferromagnet, which is simpler than the quenched one and it can 
be studied with convexity techniques which generalize those illustrated in 
[21 llOj . This study is reported in an appendix. Despite its simplicity, the 
annealed model exhibits a rich behavior. We then studied rigorously the 
quenched model using quite simple techniques, revealing interesting connec- 
tions with glassy systems. Our model is therefore a long awaited one with 
intermediate difficulty and behavior between fully connected ferromagnets 
and spin glasses. We prove that there is a high temperature region with 
identically zero magnetization, delimited by a critical line which we find, 
in which we can compute the free energy exactly. We also find and prove 
an exact expression for the free energy and for the entropy at temperature 
zero. Our approach suggests possible behaviors and techniques to deal 
with generic low temperatures, on which we will report soon 4 . Lastly, 
we briefly show some self-averaging properties of the model, which are of 
general interest in thermodynamics, and in the case of the magnetization 
are used here to control the model at temperature zero. 

We focus on the Poisson version of the model, but the main results hold 
in the Bernoulli version as well (see for more details). 

2 The model and some preliminary facts 

In this section we introduce the model and the notations, and provide some 
useful formulas which are at the basis of almost all the calculations needed 
in this article. 

2.1 Definitions 

Given a set of points, the model is defined through configurations a : i 
±1, i — 1, . . . , of Ising spins. By {iu^ju, k^, lu} we will denote families 
independent random variables all uniformly distributed on 1, ... ,7V. The 
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Hamiltonian is the random function defined by 

K 

HNicr) = — cFj^aj^ 
i/=i 

where K ]s& Poisson random variable of mean aN, for some given a G IR+, 
which is called connectivity. The expectation with respect to the random 
choice of the spins and with respect to the Poisson random variable is 
denoted by E, and it is called quenched expectation. Given a non-negative 
real number /?, the function 

^Jv(a,/3) = -^Eln^exp(-/3i?jv(CT)) 

is called pressure, and —An{oi, (3)/l3 is the free energy. Given the simple 
relation between the two, we will indifferently use either one or the other. 
The sum 

^Jv(/3) = $^exp(-/3ifjv(a)) 

(7 

is the (random) partition function, and the Boltzmann-Gibbs expectation 
of an observable O : a ^ 0{a) e M is 

When it is not confusing, wo will omit the dependence of f7 on or on 
the Poisson random variable appearing in the Hamiltonian. When we omit 
the index N in the pressure we mean to have taken the thermodynamic 
limit. The main physical quantity in this model is the magnetization of a 
configuration 

1 ^ 
= AT^"^' ■ 

i=l 

A further notation is (•) = EO(-). Throughout the paper, t G [0, 1] will be 
a real interpolating parameter, and ^ab is the Kroneker function, equal to 
one, \i A = B, equal to zero otherwise. 
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A simple calculation immediately provides the following useful form for 
the pressure 

1 ^ 

^7v(a,/3) = alncosli/3+ — Eln^ + (Ti^ (jj^ tanh /?) . (1) 

Notice that the ferromagnetism implies that the pressure of the Poisson 
model is always larger than or equal to the one of the Bernoulli version. 
Notice also that not much changes in the Poisson model if one considered a 
sort of truncated distribution, in which the probability of the integer one is 
left unchanged and the only other possible value is zero. This preserves the 
main feature of the formula at the basis of the next subsection, and allows 
to compare this truncated model with the Bernoulli one, whose pressure is 
smaller. This is therefore a way to compare the Bernoulli model with the 
Poisson one (more general considerations can be found in [15 ). 

2.2 Properties of the Poisson measure 

Given a function 5 : N ^ M and a generic Poisson variable K with mean 
^, whose expectation is denoted by E, it is easy to verify that 

E[Kg{K)] = CE[g{K - 1)] (2) 

and that 

^E[g{K)]^E[giK + l)-g{K)] . (3) 

Along the same lines it is interesting to note that the second derivative 
resembles a Laplacian 

^Eg{K) = E[giK + 2) - 2g(if + 1) + g{K)] . (4) 

These formulas are used very often in the calculations we need in the 
present work. 
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2.3 The main derivatives 

A simple use of ^ yields 



9^1 AT (a, /3) tanh (3 + fl{aigajg] 



(5) 



d(3 1 + rj(crio cr jo) tanh/3 

where the index is henceforth used to denote random variables indepen- 
dent of those appearing in fi. A simple use of ([3]) yields instead 

dAN{a, (3) ^ ^^^^ + E ln[l + n(ai„ cr,„ ) tanh /3] . (6) 
da 

These two derivatives will be constantly used in the present work. 



3 Some basic properties of the model 
3.1 Equivalent formulation 

We want to show that the Hamiltonian can be written in three different 
forms, equivalent in distribution: 

K N / \ 1,N K'J 1,N 

v—1 i—1 \ z^— 1 / ij v—1 ij 

where if is a Poisson random variable of mean aN, {K^} are independent 
Poisson random variables with mean a, {K'^} are independent Poisson 
random variables with mean a/N. Let us explain what we mean by "equiv- 
alent" from the point of view of the thermodynamics of our model. We will 
do so by means of interpolation. We will henceforth use the same for the 
Gibbs measure, even when the Boltzmannfaktor is not associated with the 
original Hamiltonian of the model, but with any generic Hamiltonian (for 
instance an interpolating one). The weights defining the Gibbs measure 
are each time easily deducible from the calculations. Take 

$jvW = -Eln^cxp/3[^a,„a,„ +^/f^V,aj] , 

£7 i,j 
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where Ki is a Poisson random variable with mean taN, and {Kq} are 
independent Poisson random variables with mean (1 — t)a/N . Compute 
now the derivative with respect to t 



1,N 

ln(l + l^((TioCTjo ) tanh 0)-^y2 ^^^^ + ^{cri<^j) tanh /?) = 



i^?^ ^ Eln(l + 17Ka,J tanh/5) " ^ E 



where the equality of the two terms is due to the expectation with respect 
to io and joi which is included in the expectation E with respect to all 
the quenched random variables (recall that io,jo are independent of the 
random site indices in the i-dependent Q). 

3.2 Convexity of the pressure as a function of the con- 
nectivity 

In this subsection we want prove the following 

Proposition 1 The pressure Ajv(a,/3) is a convex function of the connec- 
tivity a. 

Here it is useful to specify on which Poisson random variable the Gibbs 
measure depends on, we will do so by means of an index. By K we denote 
as usual a Poisson random variable of mean aN. 
Proof. Employing ([4]) one finds 

= Q^^H^ + f^xK^jo)tanh/3] 

^ ^jg^j^ [ l + ^A-+i(g»ogjo)tanh/3 " 
[ 1 + il/f (aioO-jo) tanh/3 

where K is the usual Poisson random variable with mean aN. It is also 
easy to check that 

(o-^o ajg ) + fix [cTio (Jja o'fco f^/o ) tanh /3 



1 + fix (crfc(,cr;J tanh /? 
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Prom these last two expressions, after a few calculations, one obtains 

l + QK^l(fT,„fT;o)tailll/^ 



1 + Ok (cri„ CTjo ) tanli/3 
1 + (tanh^ (3) 



2 Q^ ^i'^io^jo^ko'^lo) - ^('^io<^io)^(o'fcoO-io) 



(1 + (tanh/3)f](aioCTjo))(l + (tanh/3)f7(afe„a( J) 
Hence 

due to the ferromagnetic nature of the interactions, which implies 

and completes the proof. □ 

3.3 The infinite connectivity limit 

Recall that the Hamiltonian of the Curie- Weiss (CW) model is 

//(f>^)(a) = iW(a), 

and the associated pressure will be denoted by A^^^\p). Given n config- 
urations (T^^\ . . . , a^"\ let us define the multi-overlaps by 

N 



1 \ ^ 1 (n) 

Notice that gi = m is the magnetization. Let us prove the next 



Proposition 2 If we let a ^ oo, (3 with 2a tanh (3 = l3' kept constant, 
the pressure A^ia, (3) tends to the one of the CW model: 

lim AN{a,(3)=A'^^'^\(3') 



lot tanh ^ 



uniformly in the size of the system. 
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A more general theorem has been proven long time ago in [T]. 
Proof. Consider the following interpolation 

1 [ 11 

ijv(t) = -Eln^exp /3 ^ a,„(T,,^ + (1 - t)/?'- W (7) 

(J u—l 

between the Curie- Weiss model and its dilute version, Ki being a Pois- 
son random variable with mean taN. Clearly iAr(O) = a'-^^\p') and 
^Ar(l) = ^Ar(a,/9). A series expansion of ^ offers 

= a In cosh /3 - a J] ^(g^) tanh" /3 - ^/^'(m^) , (8) 
at — ' n 2 

n 

where the average (• • • ) depends of t because of the interpolation. When 
a oo, /? ^ with 2atanh/3 = /?' all the atanh"/? — + for n > 1, and 
the right hand side of ([5]) reduces to a tanh /3(rn^) — 13' [m?) /2 = 0. □ 
Since at the order n = \ the quenched and annealed model coincide, as de- 
ducible from the previous subsection (see the appendix) , the same identical 
proof holds for the annealed model as well, where ^{ai^ajg) is replaced by 
f](m2). 

Remark 1 The dilute model, he it quenched or annealed, reduces to the 
fully connected one in the infinite connectivity limit uniformly in the size 
of the system. 

3.4 The ground state 

From ([5]) and 

d(i N ' 

{Hn) 
ma — -— — —a 



we immediately get 
for all N, and clearly 



hm = 



The same ground state is well reproduced by the annealed model (see the 
appendix), since these formulas stay the same. 
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4 High temperature, low connectivity, sym- 
metric region 

In this section we show that at least in a region with smaU enough con- 
nectivity we can compute the free energy in the thermodynamic hmit. We 
do not prove here the existence of any symmetry breaking nor that the 
considered region spans the whole symmetric phase. These aspects will be 
studied later on, here we are about to prove only the next 

Theorem 1 In the region defined by 

2atanh/3 < 1 

the thermodynamic limit of the pressure is given by 

A{a, [3) = In 2 + a In cosh (3 . 

In particular, when a < 1/2 any value of the inverse temperature [3 fulfills 
the condition. 

Proof. First of all notice that © implies 



since the ferromagnetism imposes Q.(aif^Uj^y) > 0. As a consequence, we get 
immediately a first inequality 

AM{oi, (3) > ln2 + a In cosh /3 . 

Now consider again the interpolation ([7]) and use ([T|) to observe that 



5Ajv(a,/3) 



> In cosh (3 



da 



dAN{t) 
dt 




= a In cosh /3 
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since ln(l + x) < x and 2atanh/3 = /?'. Therefore, by the fundamental 
theorem of calculus 

AN{a,f3) < alncosh/3 + A^^^(/3') . (9) 

Now, it is well known [HI [TU] that 

lim A'-^'^\f3')=ln2 1 

and therefore putting together the two opposite inequalities we obtain 

2Q!tanh/3<l lim AN{a, f3) = ln2 + alncoshf] 

and the proof is complete. □ 

Remark 2 We obtained en passant the existence of the thermodynamic 
limit of the free energy per spin at least in the considered region. 

Remark 3 The pressure of the Viana-Bray model at high temperature is 
the same, such an expression being of quite general validity. 

This low connectivity behavior is well described by the annealed ap- 
proximation, discussed in the appendix. 

5 The model at temperature zero 

We have just seen that in a high temperature region the free energy can be 
easily computed explicitly. We will now study the model at temperature 
zero. This will also be used to prove that the described high temperature 
region cannot be extended, in the sense that outside such a region the 
solution we obtained above does not hold. 

Let us introduce the following notation for the difference between the 
pressure and its high temperature expression: 

AN{a:,(3) = ylAr(a, /3) — ln2 — aln cosh/3 

1 ''"^ 
= — Eln^]^(l + cr,crj tanh/3)^" -ln2 . 
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We have just seen in the previous subsection that in the thermodynamic 
Hmit AN{a,/3) — > if the temperature is such that 2ata.nh.f3 < 1. We will 
see in the next section that this quantity is different from zero outside the 
region. In this section we want to study the model at temperature zero, 
and we will need 

1,JV 

ij^^(a) =liminfi;v (a, /3) = -Eln V 17(1 + ^''^j)'^" -ln2 

p— >c!o iV — ' 

Let us show that the thermodynamic limit of this quantity exists. This is 
guaranteed by the next 

Lemma 1 The function NA^^\a) is sub-additive in the size of the system 
N. 

Proof. The proof can be obtained through interpolation, and it is guided 
by the reasonings described in the next subsection. Decompose the system 
into two subsystems of sizes Ni and N2, and denote by a' the spins in the 
first block, by a" the spins of the second block, while a will still denote 
configurations of the whole system. Define, for t G [0, 1], 

l.N 1,N 1,N 

m =^hiJ2 11(1 + a.a,fo + ala^fr "["[(i + ^^^^^k^ 

^ hj ij ij 

where {K^}, {Kf}, {K2}, are families of independent Poisson random 
variables with mean ta/N, (1 — t)a/Ni, (1 — t)a/N2 respectively. Denote 
by mi and m2 the magnetizations of the two blocks. A direct calculation 
gives 



< 



by convexity (the next subsection explains the mechanism at the basis of 
this result). This allows a comparison between the values of ip{t) at zero 
and at one, which describe the dependence of A^^^ on the volume involved 
in the sub-additivity. So the proof is complete. □ 

The convexity just seen will let us compute the entropy and the free 
energy explicitly. 
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5.1 Free energy 

Consider again 

1,N 

A^°\a) = ±ElnJ2l[{l + a,a,r' -ln2 . 

It is obvious that if some K^^ > 1 then only the configurations with at = 
cjj contribute. In fact the only alternative would be ct^ — —aj implying 
(1 + aiCTj)-^'^ = 0. Therefore for a given realization of the {K''^} the set 
of spins decomposes into, say, L non-empty clusters of sizes A'^i, . . . ,Nl, 
such that all the spins in a given cluster share the sign and are connected 
by non-zero links K'^^ . The quantity ^l{aiaj) can thus take only the values 
zero and one, according to whether the sites i and j are connected or not 
(i.e. belong to the same cluster or not or equivalently K'^^ is different or 
equal to zero). In the fully connected Curie- Weiss model there is only one 
cluster coinciding with the whole system of N spins, all aligned. The fact 
that O(ijja-j) is either zero or one implies 

and 

We can also write 
and notice 

V 2 j \ s^^^. , if K'i > . 

Hence all strictly positive values of the Poisson variables yield the same 
identical contribution, and therefore at least at temperature zero it is trivial 
to see that our Poisson model is equivalent to a Bernoulli one, in which the 
couplings obey K^^ = with probability po = ex.p{—a/N), and K'^^ = 1 
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with probability pi = 1 — po ~ '^k>aPk- The notation we just used 
distinguishes the BernouUi case from the Poisson one by means of the 
bar, both for the weights Pk,Pk and for the random variables K'^^,K'^^; k 
is clearly a natural number, and pk = exp{—a/N)a^ /{N'^kl). For large 
N only the dominant terms contribute, and we could equivalently take 
pQ — I — a/N, pi = a/N. Summarizing: 

Remark 4 The results we obtain regarding the entropy at temperature zero 
hold both in the case of Bernoulli dilution and in the Poisson one. For 
statements of wider validity see U5\/ . 

Notice that at any temperature the ferromagnetism implies that the Pois- 
son model gives an upper bound for the Bernoulli one, since larger values 
of K^^ (possible only in the Poisson case) increase the pressure. The oppo- 
site bound proving the equivalence between the two versions of the model is 
more involved, and relies on the fact that the mean of the coupling variables 
K^^ is proportional to so that values larger than one tend to be negli- 
gible events and the remaining two possibilities give the same contribution 
in the Poisson and Bernoulli cases. 

Remark 5 The decomposition into clusters of the spins does not depend 
on the temperature, being determined by the random couplings only. 

Moreover, it turns out that the clusters into which the systems decom- 
poses are dominated by a very large one, surrounded by many small ones, 
and this is connected with infinite dimensional percolation. 

Notice that the magnetization of the l-th cluster is Ni/N. 
The main purpose of this subsection is to prove the next 

Theorem 2 At temperature zero, the pressure per spin of the dilute mean 
field ferromagnet, be it Poisson or Bernoulli, is given by the formula 

lim lim ^ElnZAr(/3) = max{2aAf + 2aexp{-2aM) ~ aM'^} ln2 
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in the thermodynamic limit, where a is the degree of connectivity of the 
system. 

The value of M where the maximum is attained is such that 



and exhibits a critical value a — 1/2, below which M is equal to zero, 
above it is different from zero. We will get back to this along the proof. 

The statement and proof of this theorem provide a connection between 
statistical mechanics and graph theory [T^ . 

Proof. The theorem will be proven through two opposite bounds. The 
convexity we found at the basis of the sub-additivity of iV^^-* (a) makes it 
possible to introduce a "replica symmetric cavity" , whose bound will turn 
out to be exact. 

First bound. Define, for t £ [0, 1] 



where Kl^ are independent Poisson random variables with mean ta/N, Kq 
are independent Poisson random variables with mean (1 — t)2Ma with M 
a free parameter. If Kq = 0, then (1 + ai)^" = 1. If instead Kq > 0, then 
(1 + (Ti)^o = (5icr;2^o. Which means the spin cr^ is forced to take the value 
one. Hence the system decomposes into blocks with non-zero internal links 
Kl'' , and within each of these blocks all the spins are equal to one if Kq = 
all values of i, while the spins are free to fluctuate if Kq > for at least one 
value of i. Notice that n{(7iaj) = 0, 1 and il{cri) — 0, 1, moreover for t = 
the function tpN can be computed explicitly. Let us calculate its derivative 



M = 1 - exp(-2Q;M) , 



1,^ N 



^M{t) = T^Eln^ [](! + ^■''^^■)''^' 11(1 + a.)^o 



(7 ij i—1 




aEln(l + ritWz]) - 2aMEln(l + QtWi]) 
aln2(EOt[m2] - 2ME0tH) 
aln2(E0t[(m - Mf]) - aM^ In 2 



(10) 
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with an obvious meaning of Ctf Hence 

d 



^^fNit) > -aM^ln.2 . 
Integrating between zero and one 

(^jv(l) = -^ElnZjv > (i?iv(0) -aM2ln2 . 
The computation of (/5jv(0) is not difHcult, if we notice that 
5^(1 + c7i)^o = 2 if Ki = Q 
^(1 + cTi)^o = 2^0 if Ki>Q 

Therefore 

oo 

<^iv(0) = Eln^(l+cri)^« =poln2+^pfefcln2 = [exp(-2aM)+2aM] ln2 

(T fc=l 

where Pa = cxp(— 2q:M) and in the sum over the integer k the term corre- 
sponding to = can be added. Hence 

Um ^ElnZjv(/3) > {2aM + exp(-2aM) - aM^] ln2 = [<p{M)] In 2 , 

which proves the first bound for any size N of the system. Notice that the 
derivative with respect to M of <f{M) is 

<fi'{M) = 2a- 2aexp(-2aM) - 2aM = 2aM J^l^^^PhM^l _ i 

The function (1— exp(— 2aM))/M is decreasing in M and so is d(p{M) / dM^ , 
and (f>{M) is a concave function of M^. Its maximum is at zero if a < 1/2, 
it is different from zero if a > 1/2, and more precisely where 

1 - M = exp(-2aM) . (11) 

Let us recap what we proved 

Hm Hm Ajv(a,/3) > sup{2aM + exp(-2aM) - QM^}ln2 . 

/3— >oo AT— »(x) M 
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We are now going to show that the opposite bound holds in the thermody- 
namic Hmit. 

Second bound. We will make use of the self-averaging of the magnetiza- 
tion, proven in subsection 17.11 Let start from the sum rule 

lEln^(l + a.a,r^ = 

a 

[2aM + exp(-2aAf) - aM^] In 2 -f a In 2 ( EVlt[{m - M^)]dt 

Jo 

which is a consequence of (fTU)) and the expression of (pjsi (0) we computed. 
We want to show that 

lim / Ent[{m ~ M)^]dt ^ . 

Given a statement /, let us define the truth function xi as equal to one if 
/ is true, equal to zero if / is false. Let us then split the integral into three 
pieces: 

1 .1 
Ent[{m-Mf]dt = / Ent[{m- Mfxm<M~e]dt 
Jo 

+ / EQtiim - M)'^XM-e<m<M+e\dt 

Jo 

+ f Ent[{m- Mfx„,>M+e]dt . 
Jo 

The second of the three terms in the right hand side is clearly bounded by 
e^. We want to show that the other two terms vanish in the thermodynamic 
limit for any e > 0, so to have 

lim ANia,oo) <[2aM + e-xp{-2aM)-aM^]lii2 + e'^alii2 V e>0. 

We will show that the integrand vanishes for any given value of t. Consider 
the first term. In this case one has 

(to - Mf <{! + Mf 

and thus 

m<M-e\ ■ (12) 
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Similarly 

[Xm>M+e] (13) 

for the third term. Let us proceed with the first case by dividing the interval 
[— 1, M — e] into L small sub-intervals [nia, ma+i], labeled by a = 1, . . . , L. 
We assume mi = —1, m^+i — M — e. We can write 

if we define 

l,Af N 

= E 11(1 +^^^^-)''^' 11(1 + ^')''° ' 

cr ij i—1 
1,N N 

a ij i—1 

assuming ma < m < ma+i and M is chosen to satisfy pT|) . We know from 
the sum rule that 

lim AN{a,P) > [2aM + exp(-2aM) -aAf^]ln2 . 

If we knew that, choosing L sufficiently large, we also have 

^ElnZ^^j < [2aM + exp(-2aM) - aM^] In 2 



then we would be sure that 



lim = 



because of the almost certain convergence. This would mean that the right 
hand side of (fT2|) vanishes, and analogously for (fT3|) . A simple interpolation 
argument brings 

-^E In J < -^E In - taAl'^ In 2 + ta{M - rriaf In 2 , 

since r2j^[(m — Af )^] < (Al — ma)'^ if < m < ma+i < M, with an obvious 
meaning of flf. Now define 

pin) = (exp(-2Q:M) + 2aM) ln2 + (1 - M)(lncosh A - Atanh A) 
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with 



tanh A = 



H-M 
1-M ' 
M-fi 

1-M ' 



if M < fi<l , 

if 2M -l<^i<M . 



Notice that p{ii) — !■ —oo if /i < 2M — 1, and tanh A = if = M and thus 
p{M) = (exp(-2aM) + 2aM) In 2. We need at this point 

Lemma 2 Given an interval [61,62] we have 



The necessity to introduce the function p and the proof of the lemma are 
the result of a standard micro-canonical analysis. Proof. Notice that p is 
symmetrical in the interval [— 2M + 1, 1] with respect to the central point 
(1 — M, and it is increasing for 2M — 1 < p, < M, decreasing for M < p < 1. 
We are interested in computing 



1 i 
Jim — Elny'TT(l + cri)^*X6,<„<62 




1 

lim -ElnV rr[(l+ai)'^'exp(Aai)] =Elny[(l + <Ti)^'exp(A<Ti)] . 




This is easy to do as (1 + (Ti)^^ exp(A0-i) is equal to 2 cosh A if = 0, 
to 2^' exp XifK^ > 0. Hence 



Eln^[(l+ai)^'exp(Aai)] 



00 



po (In 2 + In cosh A) + pfe (fc In 2 + A) = 

k=l 

Po(ln 2 + In cosh A) + 2aM In 2 + A(l - po) 
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where clearly po = exp{—2aM) and J2k>oPk^ ~ 2aM. Let us now com- 
pute 

1 ^ 

-Eln^I](l+aO^'x„<^< 

cr i=l 
1 ^ 

^Eln^ + Gif' exp(AaO] exp(AiV(M - m)) = 

<T i=l 

A/x + exp(-2aM)(ln2 + In cosh A) + 2aMln2 - A(l - exp(-2aM)) 

since Xm<^i < exp(AA^(/x — m)) for all A > 0. But now it is easy to find the 
minimum with respect to A of this expression, which is precisely the convex 
function p previously defined: the condition is tanhA = (M — /i)/(l — M), 
as long as 2Af — 1 < /i < M , otherwise the p is decreasing and the infimum 
is for A oo. The case /U > M is analogous. This proves the lemma. □ 
Actually, standard micro-canonical approach would allow to prove the 
equality in the statement of the lemma. The lemma just proven implies 

-^E In J < p{ma+i) + ta{M - niaf In 2 - taM"^ In 2 

^ Xma<m<ma+i ^ Xm<m„+i trivially. We can rewrite the previous inequal- 
ity as 

;^IElnZ^_,< 

p{ma+i) + aln2{M -maft + aln2[{M - ruaf - (M - ma+if]t - a In 2M'^t 
< (exp(-2aM) + 2aM - aMH) + {1 - M) (In cosh A - AtanhA)-^ 

aln2(M - rrio+i)^ + 2(mo+i - rua) 

since (M - m.a)'^ - {M - ma+i)^ < 2{ma+i - nia). As tanhA = (M - 
ma+i)/(l — M), if we could now prove 

(1 - M)(lncoshA - AtanhA) +aln2(M - m^+i)^ = 

(1 - M) [In cosh A - A tanh A + a In 2(1 - M) tanh^ A] (14) 
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to be strictly negative for all ma+i < Af — e, we would be done, for it would 
suffice to take the partition {ma} fine enough in order to have 

lim In , < exp(-2aAf) + 2aM - aU'^t 

which proves the theorem. So let us consider the function 

g{X) = In cosh A - A In 2(1 - A/) tanh^ A . 

Such a function is easily proven to be concave decreasing in A^, and this 
follows from the fact that if Af (a) fulfills (dH), then it is increasing in a. 
The explicit form of the derivative of M(a) and these last observations give 



dgiX) 



dA2 

and from the concavity 



<-i(l-ln2) 

A=0 ^ 



ff(A) < -i(l-ln2)A2 < -^(1 - ln2)tanh^ A . 

Recall that tanhA = (Af — ma+i)/{l — Af), and we can now state that 
the quantity in (|14p is strictly negative whenever TOq+i < Af — e. As 
anticipated, this proves the theorem. In fact, the case to > Af + e is totally 
analogous. □ 

5.2 Entropy 

Define the entropy per spin of the model as 

SNia,P) = AN{a,P) - (idpANia, (3) . 
Its value at temperature zero is 

Sn\'^) = i™ SAr(a,/3) . 

Recall that, according to our definitions 

^Ar(a,/?) = ANia,l3) + In 2 + a In cosh /3 
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and notice 

dpAN{a,p) = a(l - tanh^/3)E < - tanh^ /3) . 

1 + S 2 ((7i„ crj(, ) tanh p 2 

Since lini^^oo /3(l—tanh^ /3) = 0, and moreover lim^^oo (In cosh /3—/3 tanh /3) 
— In2, the inequahty above means that 



We also know that 



Lemma 3 The function Ns^^\a) of the size N of the system is sub- 
additive. 



This is a consequence of Lemma [1] and guarantees that UmAr^oo Stv^Iq^) = 
inf^ s^)(a)=s(0)(a). 

Now, in the previous subsection we computed (a) in the thermo- 
dynamic Umit, and thus we also gained the next 

Theorem 3 The following formula 

s(o)(a) = sup{-a(l - M f + cxp(-2aM)} In 2 

M 

provides the entropy per spin of the model at temperature zero, both in the 
Poisson and in the Bernoulli cases. 

As expected this formula prescribes the correct critical value a = 1/2, 
below which the entropy decreases linearly from In 2 to half this value: 
s('')(q;) — {1 — a) In 2. Notice that s*^°^(q;) > for all values of a and tends 
to zero as a ^ oo. 

Somewhat surprisingly, the same low connectivity behavior with strictly 
positive entropy is captured by the annealed approximation as well (see the 
appendix). 

As a last remark, we already noticed that the free energy of the Viana- 
Bray model at high temperature is the same as the one of our model, as 
a consequence the Viana-Bray model too has strictly positive entropy at 
temperature zero. 
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6 Existence of the transition and the critical 
Hne 



Consider again 



A]\r{a,^) = Ajv(a,/3) — ln2 — alncosh/? . 



Wc saw already that in the thermodynamic Hmit AN{a,/3) 



if the 



temperature is such that 2Q!tanh/3 < 1. We want to show here that this 
is precisely the high temperature region where the magnetization is iden- 
tically zero and that in the complementary low temperature region the 
symmetric solution described above does not hold. In fact the model ex- 
hibits a transition and the magnetization fluctuates in the low temperature 
region. More precisely we want to prove the following 

Theorem 4 In the low temperature region defined by 

2a tank P > 1 

the limiting free energy differs from the symmetric expression, i.e. 

liminf Ajv(a,/3) > . 



N—^fX 



Proof. Let us keep 2atanh/3 = f3' constant, so that 



2(tanh/J)da + 2a(l - tanh^ /3)d/3 = . 



We have 



dAN{a,P) = daAN{a, P)da + d0A]\[{a, P)d(3 
= Eln(l + n{ai„(Tj„) tanh/?)rfa 



(15) 




dl3 



= E ln(l + 0((7io aj„ ) tanh /?) - E 



da 



1 -I- ri(crj(,crjo)tanh/? 
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which is non-negative as ln(l + x) > x/{l + x) for x > 0, and we are 
considering the case da > 0. The function hi(l + x) — x/{l + x) is positive, 
increasing and convex for < a; < 1. This allows for the computation of 
Ais[{a, (3), at least for — > oo and /3 ^ 0. 

As anticipated, our strategy for the proof consists in considering the 
limit as — > oo, so recall the notation 

= lim liminf i7v(a,/3) = — ElnV TT(1 + ct^CTj)^" -ln2 

where once again K^^ are independent Poisson random variables with mean 
a/N, and we know that 



hm i^^ (a) = for a<- 



Lemma [T] implies limAr^oo ^nH'^) — ii^^Af ^n\'^) = A^°\a), and Theorem 
M gives 

A^°\a) = nmx{2Ma-|-exp(-2aM) -Q:M2}ln2 . 



Notice that 

-^^^^'^""^ 

which means A'^*'^(a) > for a > 1/2, and finally implies the statement of 
the theorem we wanted to prove because of ([T5|) . □ 

Remark 6 In this last lemma the maximum is clearly where and only 
where 1 — M{a) = exp(— 2Q!M(a)), which determines a — 1/2 as the crit- 
ical point. In fact, M{a) = for a < 1/2, while M{a) > for a > 1/2. 
Furthermore, the critical index is given by M{a) = 2y/a-l/2 + 0{M^) 
for a > 1/2. 

Notice that the annealed model has the same high temperature free 
energy and the same critical line as the quenched one (see the appendix). 

The interested reader might want to compare our proof with the results 
obtained through replica techniques in [13] . 
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7 Self-averaging 



In this section we briefly discuss the hmiting self-averaging of the free 
energy density, of the magnetization, and we exhibit the relations coming 
from the self-averaging of the internal energy. The self-averaging properties 
hold for both the quenched and the annealed models (in the latter the 
relations simplify as the squared magnetization replaces the product of two 
randomly chosen spins). 

7.1 Magnetization 

The magnetization turns out to be self-averaging, and therefore it is an 
actual order parameter, although it is not clear whether it is the only order 
parameter. Let us prove the next 

Theorem 5 The following identity 
holds in the thermodynamic limit. 

Proof. Consider the derivative with respect to the connectivity of the 
mean squared magnetization. Using ^ we have 

—En(m^) = NE\nK+i(m^) - f^^fm^)] 
da 



which is bounded. Proceeding further 
da 1 -I- f](crioCrj„)tanh/3 



iV(tanh/3)E- ^(-''^-^^o) - ^im^Ma^^a^o) 



(1 + A(o)(a)17(cr^oCrjo) tanh/3)(l-f r2(m2) tanh/3) 
+Ar(tanh' ^)]g_^^("^')[^^("^''^»o'^jo) - ^im^M<y,„aj,)] 



(1 + n{aia(Jja)tanhp){l + n{m'^) tunh l3) 



but this means 

E[fl{m^a,^aj^) - n{m'^)n{a.,„aj„)] = En{m*) - En^{m'^) 
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when A'^ — oo. □ 

Another even stronger self-averaging identity is formulated in the next 
Theorem 6 The following identity 

holds in the thermodynamic limit. 

We are not going to prove this here, we will show some similar relations 
later in this section. 

7.2 Free energy 

As expected on a physical ground, the free energy of our model does not 
fluctuate in the thermodynamic limit. This is formalized in the next theo- 
rem. 

Theorem 7 For all values of N , a, and (3, the following exponential in- 
equality 

> e| < 2 exp[7V(e - a(l -f e/a) ln(l + e/a))] 

estimates the probability for the random free energy to deviate from its 
expectation. 

Sketched proof. The proof is quite standard and proceeds along the 
same lines dictated for instance in [16J for disordered systems. We will 
only sketch the proof here. For a given real number A, let us define 

= In El exp[AE2 InZjv(t)] 

with 

cr ^ l/=l l/=l 

where Ki and K2 are Poisson random variables with mean taN and (1 — 
t)aN respectively, while Ei and E2 denote the expectation with respect to 



In Zm E In Zm 

I3N I3N 
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{Ki, {ii/}, {ju}} and to {K2, {fci^}, {1^}} respectively. The reason for such 
a construction is 

exp[$^(l)-$^^(0)] =Eexp[A(lnZAr-ElnZjv)] . 

The derivative with respect to t reads 



a at 



l.N 

E 



El exp[AE2 lnZ7v(t)(exp(AE2 In f](cr,(Tj) - 1) - AE2 In ^(crifTj))] 



Eiexp[AE2 InZAr(t)] 
The simple observation that 

-P < E2lnr2[exp(/3cricrj )] < P 



ensures 



< aiV[exp(|A|/5) - 1 - |A|^] 



The validity of this bound for all A together with Tchebyshev's inequality 
implies that the random free energy deviates from its expectations with a 
probability exponentially small in the size of the system. □ 



7.3 Internal energy and self-averaging relations 

The self-averaging of the free energy implies, by standard convexity argu- 
ments [5], the self-averaging of the internal energy. More explicitly, let us 
state without proof the next 

Theorem 8 In the thermodynamic limit the internal energy does not fluc- 
tuate 



lim 



N 



N 







with the possible exception of a zero measure set of values of the inverse 
temperature [3. 
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From the self-averaging of the internal energy one can obtain useful infor- 
mation of probabilistic nature about the thermodynamic behavior of the 
model. Since we already have an expression for the averaged internal en- 
ergy, given by ([5]), it is not difficult to perform calculations similar to all 
the others done so far and obtain 

Proposition 3 In the thermodynamic limit, the following identity holds 

^ fl{(7i(7jak<Ji) + Q{ai(7j) tanh/3 + fl(ak<Ji) tanh f3 + tanh^ /3 
[l -I- fl{aiaj) tanh/? -I- fl{akO'i) tanh/3 + Q{ai(7jakO'i) tanh^ /3_ 

_ tanh/3 -I- r2((7iCTj) ^ 
1 -t- ri((Tfe(T;) tanh/3 _ 

for all i, j, k, I £ {1, . . . , N}. 

This identity provides a further restriction beyond the self-averaging of the 
magnetization, characterizing the thermodynamics of the model. 

8 Outlook 

We showed that dilute mean field ferromagnets are not so poorer than 
spin glasses, but not as difficult either, and this makes them quite inter- 
esting. The annealed model is already rather interesting, contrarily to 
other annealed models. It has a phase transition, and a non-negative zero- 
temperature entropy. The annealed model actually enjoys a strictly posi- 
tive entropy at temperature zero. It provides the same ground state as the 
quenched model, and the same high temperature regime, with the same 
critical point. The quenched model is highly non-trivial. The control of 
the zero-temperature regime we gained suggests an approach to the model 
at a generic temperature, since the kind of "replica symmetric" method 
employed here allows for an extension to a more general distribution of the 
order parameter fl{aiaj). The approach based on cavity fields developed 
instead in the physical literature for spin glasses, though not fully rigorous, 
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also can be used to study the generic behavior. We will report our results 
on these topics elsewhere [4], along with other interesting developments: a 
(G)REM-like approach, the p-spin case, the form of the cavity fields and its 
consequences. Both the annealed and the quenched models enjoys the sta- 
bility properties known for spin glasses: a suitable cavity field is equivalent 
to the addition of one particle to the system, and this makes all overlaps 
squared. As a consequence, the free energy is the difference between two 
terms, and does not depend on certain overlap monomials, and this makes 
it possible to compute the critical exponents. This is all quite easy in the 
annealed model, while in the quenched model we have no proof, through 
our techniques at least, of the existence of the thermodynamic limit of the 
free energy density as yet, and the procedure to compute the critical expo- 
nents in slightly more involved anyway. We will report on this in a separate 
paper [1] . As a last remark, some generalizations to a bipartite model with 
two populations of spins are possible [3]. 

A Annealed model 

In this appendix, we have to consider separately the expectation with re- 
spect to the Poisson variables and the expectation with respect to the 
random choice of the spins. The expectation with respect to the random 
choice of the spins is denoted by Es , while the expectation with respect to 
Poisson random variables is denoted by Ep. Collectively we will use again 
E — KpKg. The annealed pressure is defined by 

AN{a,(3) = lEplnE,^cxp(-/?i/Ar(a)) . 

cr 

We may easily include the action of an external field on the system, by 
adding a term —hNm to the Hamiltonian. The next two subsections are 
devoted to the calculation of the limiting pressure by means of two opposite 
bounds. Both bounds are obtained following the ideas of [HITO], which 
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allow for a generalization when a quadratic function is replaced by a generic 
convex one. 

It is easy to see that the pressure can be written as 
AN{a,l3) = a In cosh /3 + -^Ep In^ exp(inn(l + tanh/3)) + N/3hm). 

It is thus convenient to define a function / of the magnetization m by 

/(m) = ln(l + n? tanh /3) , (16) 

whose derivative is 

777 

/M^2tanh/3 ^^^,^^^^^ (17) 
and its expression will be used often in the rest. 

A.l Lower bound for the annealed pressure 

Notice that the function / defined in is easily verified to be convex, 
and therefore 

/(to) >/(M) + /'(M)(77i-M) 

for any given M. This means that we can proceed like in [9l[T0], where the 
function was /(to) — rr? ^ and prove the next 

Lemma 4 The following bound 



M 

(2a tanh B) ^ 



Ajv(a,/?) > sup<^ ln2 + alncosh/? + aln(l + Af'^tanh/?) 

' (2a tanh B) hln cosh 

^ ^ 'l + A'P tanh (3 



holds for all values of the size N of the system. 

Proof. From the just mentioned convexity of / it is obvious that 

A]y{a,f3) > a In cosh /3 

+ sup I ^Ep \n^cxp{K[f{M) + f'{M){m - M)] + PNhm) 



M 
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This is the lower bound. We will proceed a few steps forward to have a 
more explicit expression of the annealed pressure. Let us define 

Atriai{M) = aln cosh/3 + -^Ep In^ exp(is:[/(Af) + f'{M){m - M)] + /3iVm) 

= alncosh/3+ ^Epln{exp(i4:[/(M) -/'(M)]) 
x'^exp[{Kf{M)+ph)Nm]} 

= ln2 + alncosh/3 + a[/(Af) -/'(Af)] 
+Ep In cosh 



^f'(M)+ph 



Using again a simple convexity argument we take the expectation Ep inside 
the In cosh function and obtain 

Atrtai{M) > ln2 + alncosh/3 + a[f{M) - /'(Af)] + lncosh[a/'(Af) + /3/i] 

although the equal sign would hold in the thermodynamic limit since K/N —> 
a. Using the explicit form of the function / we can write 



AtriaiiM) > ln2 + alncosh/3 + Q;ln(l + A'f^tanh/3) 



-(2Qtanh/3)- 



1 + A/2 tanh (3 
The final result is thus 



~ln cosh 



(2a tanh /3) 



M 



1 + Af 2 tanh 13 



I3h 



AN{a,P) > supi ln2 + alncosh/3 + aln(l + Af^tanh/?) 
M L 

Af2 

- (2a tanh /3) z rr?; -. — - + In cosh 



1 + A.f2 tanh /3 



(2a tanh /3) 



M 



1 + A/2 tanh /3 



for any size N of the system, which is precisely the statement we wanted 
to prove. □ 

A. 2 Upper bound for the annealed pressure 
Like in [HI [TD] , we will employ the following trivial identity 

SmM = 1 



M 
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to prove the next 



Lemma 5 The following bound 



AN{a,l3) < ^'^^^^^ + sup I ln2 + alncosh/g + Q;ln(l + tanh^) 



M 



M2 

- (2a tanh B) h In cosh 



M 



holds for all value of the size N of the system. 
Proof. One has 

^7v(Q!,/3) = alncosh/3+ -^Ephi^^5„Mexp[iirin(l + TO^tanh/3) 

M a 

+Nphm] 

= alncosh/3+ -^Epln^^(5„A/exp[ii:/(TO)] +iV/3/iTO] , 

M a 

where / is again the one defined in (fTH)) . But now thanks to the deha 
function 

f{m) = f{M) + f\m - M) 

so that 

AN{a, (3) = a In cosh /3 

+ ^Ep hi^ ^ 5^M e^v[K{f{M) + f'{M){m - M))] + N (5hm] . 

M <y 

At this point we trivially observe that 

and thus 

A^ia, /?) < alncosh/3 

+ -^Ep In^ ^ cxp[is:(/(Af ) + f'{M){m - M)) + N f3hm] . 

M a 
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Observe now that M can only take 2A'' + 1 values, therefore 



AN{a,P) < a In cosh /? 

+ sup ^mn[{2N + 1) y exp[i('(/(M) + /'(m - M)) + Nphm] 
M N ^ 

which means 

- , ln(2iV+l) 

^iv(Q!,/?) < -^-^j ^ + SUpAtria!(M) , 

iV M 

according to the definition of An-iaiiJ^l) given in the previous subsection. 
Therefore we have proven the lemma. □. 

A. 3 The annealed pressure 

We can put together the lemmas of the previous two subsections, and 
summarize the final result in the thermodynamic limit. 

Theorem 9 The limiting annealed pressure is given by the formula 



.1 



lim AN{a, P) = sup< ln2 + alncosh/3 + aln(H- M^tanh/3) 



M2 

- (2a tanh (3) —5- — +ln cosh 

^ ^l + M^tanh/? 

for all values of a, (3 and h. 



M 

(2a tanh/?)- — r^+Ph 



As we said, this follows immediately from the lemmas proven in the previous 
two subsections, which together imply 

lim AN{a,f3) = sup Atriai{M) 

and noticing again that K/N a one has statement of the theorem. 

Remark 7 Notice that the convexity of the function f allows to prove the 
existence of the thermodynamic limit for the free energy per spin apart from 
its calculation, using now standard techniques. 
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A. 4 Symmetric region 

Let us remove the external field by taking h ^ 0. 

In the expression of the annealed pressure just computed, the supremum 
lies where 



d 
dM 



af(M) ~ af'{M)M + In cosh a/' (Af) = 

a/"(M)[tanh(a/'(M)) - M] = 

where / is once again defined in (jl6p . This implies 

M 



tanh 



(2atanh/3)- 



M 



1 + M2 tanh/3_ 

from which we deduce the critical point discriminating the region where 
there is only the zero-magnetization solution M — from the region where 
there are two non-zero opposite solutions ±M*: 

2a tanh /? = 1 

is the condition that defines the critical line. So thanks to the infinite 
connectivity limit one can deduce the critical point of the finite connectivity 
annealed model from the critical point of the fully connected model. 
When AI = the annealed pressure is 



N 



lim ylAr(a, /?)= In2 + aln cosh/3 



which holds where 2atanh(/3) < 1. Therefore the annealed and quenched 
pressures coincide in this region. 



A. 5 Positivity of the entropy at temperature zero 

Recall that 

Definition 1 The limiting entropy is 
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The value sq the entropy takes at temperature zero is 

so{a) = lim A{a,P) - l3dfjA{a,(}) . 

Notice that when (3 ^ oo the condition for the supremum in the formula 
of the annealed pressure is attained at 

M = Oifa<-; M = tanh if a > ^ . 

- 2 1 + M2 2 

When a < 1/2 we know 

A(a,/3) =ln2 + alncosh(/3) V /? 

and it is easy to sec that ln2 > SQ{a) = (1 — a) In 2 > (ln2)/2 > 
in this case. Notice the high degeneracy of the ground state: when the 
connectivity is too low there are too few interactions to move the entropy 
away from the value In 2 it has when the absence of interactions makes all 
configurations equally probable, maximizing the entropy. 

Let us assume a > 1/2. When computing a derivative of the pressure we 
do not have to differentiate it with respect to M, because of the supremum 
condition. Notice also that increases with a, and in particular 1 
as a ^ 00. Moreover, simple calculations yield 

so(a) = (1-a) ln2+aln(l+M2)-2a^^^+lncosh (^2a^^^^ (18) 

which is even in M as expected (and it is (1 — a) In 2 for M = 0, recovering 
the case of low connectivity). This expression for so{a) also says that the 
zero-temperature entropy tends to zero when a — > oo. Now 

°^ ' = -ln2 + ln(l + Af2)-2-— - + -— -tanh 



da • y ■ I 1 + M2 1 + M2 1 + M2 

= -ln2 + ln(l + M2) < 

Notice that such a derivative tends to zero when a — > oo, since in this limit 
M ^ 1, and it is increasing in a. In other words the zero- temperature 
entropy is convex in the connectivity. 
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So the zero-temperature entropy, as a function of the connectivity, de- 
creases linearly from In 2 to (ln2)/2 for < a < 1/2, then it becomes 
strictly convex, decreases monotonically, and asymptotically decays to zero. 

In particular, we proved the following 

Proposition 4 The entropy of the infinite volume annealed model remains 
strictly positive 

so{a) = lim [A{a,/3) - pd0A{a,l3)] > 
when the temperature decreases to zero. 

Remark 8 We encountered an annealed model without the problem of a 
negative entropy: on the contrary, it has a strictly positive entropy with 
a highly degenerate ground state, and the entropy vanishes only when the 
number of interactions is large enough, namely when a — > oo. 

Notice that in the low connectivity region the annealed entropy is the 
same as the one of the quenched model. 

A. 6 Infinite connectivity limit 

We already saw that both the quenched and the annealed model reduce to 
the fully connected one in the infinite connectivity limit. The same result 
can be obtained directly from the formula of the free energy we found (only 
in the thermodynamic limit) which reduces in the infinite connectivity limit 
to the well known formula for the free energy of the Curie- Weiss model. 
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